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Abstract 

\ Results for penetrating probes treated within a hybrid hydro-kinetic model are pro- 

' jected onto the energy range covered by the NICA and FAIR projects. A new source 

PsJ . of dileptons emitted from a mixed quark-hadron phase, quark-hadron bremsstrahlung, 

Q_ii is proposed. An estimate for the vrvr — a ^ 77 process in nuclear collisions is given. 

^ . 

D.I. Blokhintsev, whose centennial anniversary of the birthday this conference is de- 
voted to, has contributed to various fields of physics and its apphcations, particularly, to 
hydrodynamics. Sixteen years ago, at the dawn of hydrodynamics he made an important 
remark [1] concerning possible violation of the uncertainty principle in initial conditions of 
O ■ the Landau hydrodynamic theory. From up-to-date view, this Blokhintsev's estimate looks 
slightly naive but in principle it is correct until now and should be taken into account in 
modern development of a relativistic hydrodynamic approach. Here we present a hybrid 
model which combines kinetic and hydrodynamic descriptions. To certain extent it can be 
!>■ \ considered as a possible solution of the problem put by D.I. Blokhintsev. 

In the hybrid model [2|, the initial stage of heavy ion collisions is treated kinetically 
tJ" ■ within the transport Quark Gluon String Model (QGSM) [3j whereas the subsequent stage 
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is considered as an isentropic expansion of a formed dense and hot system (fireball). The 
transition from one stage to another is solved by considering the entropy evolution. 
00 \ In FigH], the ratio between entropy S and baryon charge Qb oi participants is shown 

^ ' for In+In collisions at the impact parameter 6 = 4 fm and bombarding energy 158 AGeV. 

Being calculated on a large 3D grid, this ratio is less sensitive to particle fluctuation as 
^ . compared to the entropy itself. Small values of the baryon charge at the very beginning 
of collision result in large values of the S/Qb ratio. It is clearly seen that for tkin ^ 1-3 fm/c 
this ratio is practically constant and this stage may be considered as isentropic expansion. 

To proceed from kinetics to hydrodynamics, we evaluate conserved components of the 
energy- momentum tensor Too, Tqi, T02, T03 and baryon density Hb (the zero component of 
the baryon current) within QGSM at the moment tkin = 1-3 fm/c in every cell on the 3D 
grid. This state is treated as an initial state for subsequent hydrodynamic evolution of a 
fireball. The time dependence of average thermodynamic quantities is presented in the left 
panel of FigO 

The latter stage is evaluated within the relativistic 3D hydrodynamics [2J. The key 
quantity is the equation of state. In this work, the mixed phase Equation of State (EoS) is 
applied [5] which allows for coexistence of hadrons and quarks/gluons. This thermodynam- 
ically consistent EoS uses the modified Zimanyi mean-field interaction for hadrons and also 
includes interaction between hadron and quark-gluon phases, which results in a crossover 
deconfinement phase transition. In addition to [5J, the hard thermal loop term was self- 
consistently added to the interaction of quarks and gluons to get the correct asymptotics 
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Figure 1: Temporal dependence of entropy S per baryon charge of participants (left 
panel) for a semi-central In+In collision at Eiab =158 AGeV. In the right panel the average 
energy (solid line) and baryon (dashed) densities of an expanding fireball formed in this 
collision. Dotted line shows a contribution of quarks and gluons to the energy density. 

at T >> Tc and reasonable agreement of the model results with lattice QCD calculations 
at finite temperature T and chemical potential [S]. This agreement is demonstrated 
in Figl2] where the reduced pressure Ap/T^ = {p^fis) — p{f^B = 0))/T'* is compared with 
recent lattice QCD data. 

The fraction of unbound quarks/gluons defined as Ppi/p = {riq + + ng)/{nq + + 
ng + ub + um) is presented for the mixed phase EoS in the right panel of Fig|2J It is seen 
that even at a moderate temperature T ~ 50 ^ 100 MeV the quark/gluon fraction sharply 
increases at the baryon density ^b/jt-o ~ 6 and dominates thereafter. At T = 200 MeV the 
admixture of hadrons is, naturally, quite small. 

Consider now penetrating probes. To find observable dilepton characteristics, one should 
integrate the emission rate over the whole time-space x = {t, x) evolution, add the contribu- 
tion from the freeze-out surface ('hadron cocktail'), and take into account the experimental 
acceptance. To simplify our task, we consider only the main channel vrvr ^ p ^ l^l^. In 
this case the dilepton emission rate is 



where the integration is carried out over the whole space grid and time from t = till the 
local freeze-out moment. Here = = qQ—q^, fsilo, T{x)) is the Bose distribution func- 
tion, and C{M) is the lepton kinematic factor. The imaginary part of the electro- magnetic 
current correlation function lmIlern{(l,T{x), fib{x) includes in-medium effects which may be 
calculated in different scenarios. The recent precise measurements of muon pairs [7] allowed 
one to discriminate two main scenarios, in particular those based on the Brown-Rho (BR) 
scaling hypothesis [8] assuming a dropping p mass and on a strong broadening of p-meson 
spectral function as found in the many-body approach by Rapp and Wambach [9]. It was 
shown [7J that the measured excess of muons is nicely described by the strong broaden- 
ing of the p-meson spectral function. In contrary, the BR scaling hypothesis predicts a 
large shift of the p-meson maximum towards lower invariant mass M in contradiction with 
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Figure 2: Temperature dependence of the reduced pressure (left panel) at the baryon 
chemical potential fiB = 210, 330, 410 and 530 MeV (from the bottom) and fraction of 
unbound quarks (right panel) within the mixed phase EoS. Points are lattice QCD data for 
the 2+1 flavor system [6]. 



experimental data [7]. 

This result looks quite disappointing. First, one sees no signal of a partial restoration 
of the chiral symmetry for the sake of which dilepton measurements were originally un- 
dertaken. Second, to be consistent with the QCD sum rules both collision broadening and 
p-mass dropping should be taken into account [lOj. Generally, the exact relation between 
the p mass and quark condensate is not flxed by the QCD sum rules in contrast with the 
BR scaling, so it is questionable whether to prescribe some T dependence to the BR mass 
shift. This shortcomings of the analysis in [7] are commented by Brown and Rho [11]. In 
this respect, for the pion annihilation tttt ^ p — > we estimated [2] the imaginary 

part of the p-meson self-energy in the one-loop approximation assuming the Hatsuda-Lee 
relation for the modifled p mass: m*(x) = mp{l — 0.15 • nB{x)/no). The calculated result is 
presented in Figl3]by the solid line |2]. Indeed, the shift of the p-meson spectral function is 
not so drastic as in [7]. In addition, one should note that dileptons carry direct information 
on the p meson spectral function only if the vector dominance is valid [11]. It is not the 
case in the Harada-Yamawaki vector manifestation of hidden local symmetry p2] . 

As seen from FigI31 the muon yield is underestimated at both low and high values of 
the invariant mass M. In the broadening scenario a low M component is explained 

by particle-hole excitations but not a partial chiral symmetry restoration. However, the 
existence of the mixed quark-hadron phase, those study is the main aim of the Nuclotron- 
based Ion Collider fAcility (NICA) project [13], may give rise to a new dilepton source, 
quark(antiquark)-hadron bremsstrahlung. Similarly to the np bremsstrahlung, the process 
for an antiquark-hadron collision may roughly be estimated in the soft-photon approxima- 
tion as 



Here the averaged cross section a{s) = al^ [s/{mN + raq)'^ — 1] and the elastic qN cross 
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Figure 3: Invariant mass distribution of dimuons (left panel) from semi-central In+In col- 
lisions at the beam energy 158 AGeV. Experimental points are from [7]. The solid line 
corresponds to the T-independent dropping mass [2] and the dashed one is the contribution 
of the quark-hadron bremsstrahlung channel. In the right panel one compares dynamical 
trajectories projected onto the T — ub plane for central In+In(158 AGeV) and Au-|-Au(40 
AGeV). The shaded region roughly corresponds to the hadronic phase. 



section is approximated by the quark scaled NN cross section [14j 
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So the production rate will be 

-m^ = J'^''J Wf^^ -^'^^^ I (2^0^^^^^'^^^)^ ^M^^^'^) ' 

where Vrei is the relative velocity of colliding qN particles and the integration in (jl]) should 
be carried out over the whole space-time available for the mixed phase. The free mass 
is used for a nucleon and uiq = 150 MeV for an antiquark. In a real case one should 
also add contributions from all other baryons, as well as that from quark-antibaryon and 
quark(antiquark)-meson interactions with proper cross sections. All these uncertainities 
are effectively introduced in ([2]) by an arbitrary factor K. 

As follows from the dashed line in FigJSl at rather reasonable value oi K = 10 the 
qN bremsstrahlung source improves agreement with experiment. It is of interest that the 
contribution of this new source decreases when the bombarding energy goes down till the 
NICA energy range (< 40 AGeV) while the contribution from particle- hole excitation is 
expected to grow since the baryon density in this range is higher (see below) allowing, 
in principle, disentangling of these two sources. The underestimated yield at high M, 
intermediate mass dileptons, can mainly be described by Drell-Yan process in the quark 
phase [15]. The mixed quark-hadron phase should also contribute to the intermediate M 
region. Its contribution can be taken into account in the way as the Drell-Yan process 
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with an additional hadron form factor. Effectively, it will increase the Drell-Yan lepton 
yield [15]. 

The phase distribution of all space cells at the early evolution moment t = 0.3 fm/c 
projected on the T — ub plane is presented in Figl31 It is seen that at the maximal NICA 
energy the baryon density in the hadronic phase for central Au+Au collisions is noticeably 
higher than that at the SPS energy in In+In collisions. It means that the difference between 
the BR scaling and broadening scenarios is expected to be more pronounced at the NICA 
energy. 




Figure 4: Quark diagrams for the 27 production in the Born approximation (a), (b) and 
through the a resonance (c). 

Annihilation of two pions into two photons is of particular interest since its cross sec- 
tion is sensitive to changes of the a-meson properties which occur in the vicinity of chiral 
restoration phase transition. With increasing temperature and density the a meson changes 
its character from a broad resonance with a large decay width into two pions to a bound 
state below the two-pion threshold m„{T,iiB) ^ 2m^(T, /i^). The calculation of the pho- 
ton pair production rate at the given T as a function of the invariant mass shows a strong 
enhancement and narrowing of the a resonance at the threshold due to chiral symmetry 
restoration [16] . We make the first estimate of this channel for a particular nuclear collision. 
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Figure 5: T-dependence of a and vr masses (left panel) and total (dashed line) and resonance 
(solid line) invariant mass distributions of two photon pairs created in the central Au+Au 
collision at 40 AGeV (right panel). 

In- medium tttt 77 process is evaluated within the NJL model [IE]. Besides the 
resonance diagram, the dominating Born terms are considered, see FigJH So the total 
rate has the Born term, resonance term and interference between them: dNiot/dM"^ = 
dNll^^/dM'^ + dN^^JdW^ + dN^Zerf/dM^ which should be substituted in eq.(^. 
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The temperature dependence of a and vr masses for this model is shown in Fig|5l The 
regime m„(T,fiB) < 2m^(T, /i^) starts at T ~ 165 MeV. The photon yield as a function of 
photon pair invariant mass is presented in FiglS] for central Au+Au collisions at 40 AGeV. 
The total number of photon pairs sharply increases above the threshold 2m.,^ and then 
flattens on the level of ~ 10~^. The resonance channel of interest is lower by about the 
order of magnitude as compared to the total yield and exhibits a spread weak maximum. 
The maximum predicted in [16] for the fixed T is washed out, as seen from Fig|5l It is not 
an easy but promising experimental problem to select out this maximum from the total 
distribution. Note that such an analysis should be carried out on a huge background of 7 
decays of 'hadron cocktail'. 
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